Focusing through random media in space and time: a transmission 

matrix approach 
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We exploit the evolution in time of the transmission matrix following pulse excitation of a random medium to 
focus radiation at a selected time delay t' and position r. The temporal profile of a focused microwave pulse 
is the same as the incident Gaussian pulse. The contrast in space at time t' of the focused wave is determined 
by the participation number of transmission eigenvalues M' and the size N' of the measured transmission 
matrix. The initial rise and subsequent decay in contrast observed reflects the distribution of decay rates of the 
quasi-normal modes within the sample. © 2013 Optical Society of America 

OCIS codes: (290.4210) Multiple scattering; (290.7050) Turbid media; (030.6140) Speckle. 



A short pulse incident upon a random medium is scram- 
bled to produce a multiply peaked speckle pattern of 
intensity in space and time. The scattering of light 
severely compromises potential applications in imaging, 
two-photon microscopy, nano-surgery and telecommu- 
nications [1-4]. The ability to focus a pulse after it is 
transmitted through a random medium was first demon- 
strated in acoustics by time reversing the transmitted 
signal [5]. The transmitted pulse is picked up by an ar- 
ray of transducers and then played back in time. The 
signal converges to a focus at the location of the source 
antenna. Recently, optical pulses have also been focused 
through random samples by shaping the wavefront [6] 
using a genetic learning algorithm based on feedback of 
light arriving at a selected point and time delay [7] or 
via nonlinear feedback of two photon fluorescence [8]. 

Transmission through a random system is fully de- 
scribed by the field transmission matrix t, whose ele- 
ments tb a are the transmission coefficients of the field 
between incoming and outgoing channels, a and b, re- 
spectively [9, 10]. Measurements of t have been used to 
focus monochromatic light [11] and microwave radia- 
tion [12] through random samples via phase conjuga- 
tion [11-13], in which the electric field from different 
points a on the incident plane arrive in phase at the 
target to interfere constructively at a maximum value 
of focused intensity. The spatial contrast for monochro- 
matic illumination is determined by the measured eigen- 
channel participation number of the transmission matrix 
M' = (J2n r «) 2 / J2n T 'n an d the size N' of the measured 
transmission matrix [14]. Here T n are the eigenvalues of 
the matrix product tv. 

In this Letter, we show that pulsed radiation can 
be focused in space and time by phase conjugating 
a time dependent microwave transmission matrix at a 
selected time t'. The temporal profile of the focused 



pulse is the square of the field correlation function in 
time. This is independent of time and for a Gaussian 
incident pulse is equal to the profile of the incident 
pulse [15]. The spatial contrast at t' is determined by 
the eigcnchannel participation number M'(t') and size 
N' of the measured transmission matrix at that time, 

M'(t') = (E„T«(i'))7E n T-nv*0, in wh ich the sum of 
T n (t') yields the transmission at that time. We find that 
as soon as a random speckle pattern forms in transmis- 
sion, the time evolution of the contrast in focusing equals 
1/(1/M'(i') — 1/N'). The contrast in focusing reaches a 
peak and then decays as a result of the time variation of 
the effective number of quasi-normal modes contributing 
to transmission. 

The samples studied were composed of alumina 
spheres with a diameter of 0.95 cm and refractive index 
3.14 embedded in low-index Styrofoam shells. The ran- 
domly positioned spheres fill in a copper tube of diameter 
7.3 cm at an alumina filling fraction of 0.068. The field 
transmission matrix t is measured for microwave propa- 
gating through random waveguides with length L = 61 
cm over the frequency range 14.7-15.7 GHz in 3200 steps. 
The number of propagating waveguide modes N of the 
empty waveguide varies from 64 to 72 over this frequency 
range. The field transmission coefficients for radiation 
polarized along the detecting antenna is measured be- 
tween all pairs of source and detection antennas posi- 
tioned on a grid of iV'=45 points on the input and output 
surfaces of the sample tube. The grid spacing of 9 mm is 
close to the field correlation length at which the field cor- 
relation function passes through zero. New statistically 
equivalent samples are produced by momentarily rotat- 
ing and vibrating the copper tube after t is recorded. 
The sample studied is diffusive with an average value of 
the steady-state transmittance (T) ~ 4. (T) is equal to 
the dimensionless conductance g. Weak absorption in the 
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alumina spheres and in reflection from the copper tube 
reduces transmission by 25%. The impact of dissipation 
upon the measured contrast is reflected in a reduction in 
M' in the steady state measurement of the order of 2% 
in the samples studied. 
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Fig. 1. (Color online) Spatiotemporal control of wave 
propagation through random waveguide, (a) Typical re- 
sposnc of h a (t') and the time of flight distribution 
{I(t')). The incident pulse is sketched in blue curve, (b) 
and (c), phase conjugation is applied to focus the mul- 
tiply peaked pulse ha(t') shown in Fig. 1(a) at t'= 33 
ns and 40 ns at the center of the output surface, re- 
spectively. The Whittaker-Shannon sampling is used to 
obtain high-resolution spatial intensity patterns shown 
in the inset of (b) and (c). 

We obtain the time dependent transmission matrix 
from spectra of the transmitted field between all points 
a and b, tb a (i>). These spectra are multiplied by a 
Gaussian pulse centered in the measured spectrum at 



vq =15.2 GHz with bandwidth a v =150 MHz and then 
Fourier transformed into the time domain. This gives 
the time response at the detector to an incident Gaus- 
sian pulse launched by a source antenna with bandwidth 
<Tf=l/27r/cr, / . The time variation Ib a {t') of an incident 
pulse launched at the center of the input surface and 
detected at the center of the output surface in a sin- 
gle realization of the random sample is shown in Fig. 
1(a). Individual peaks in intensity have widths compa- 
rable to the width of the incident pulse. The average of 
the time of flight distribution (1(f)) over transmission 
coefficients for N' 2 pairs of points and 8 random con- 
figurations is also presented and seen to be significantly 
broadened over the incident pulse. 

In order to deliver radiation to a point at the center 
of the output surface of the waveguide /3=(0,0) at a se- 
lected time t\ the incident field is constructed by phase 
conjugating the transmission matrix at time t' , t% a (t'). 
The results for t'= 33 ns and 40 ns arc shown in Figs. 
1(b) and 1(c), respectively. In both cases, a sharp pulse 
emerges at the selected time delay with intensity peaked 
at /3=(0,0). 
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Fig. 2. (Color online) Profile of the focused pulse com- 
pared with square of the field correlation functin in time 
and the profile of incident Gaussian pulse. The peak of 
the focused pulse shown in Fig. 1(b) has been shift to 
Ai=0 ns. All curves are normalized to unity at Ai=0 ns. 

The spatial profile of focusing for a monochromatic 
wave above an enhanced constant background is equal 
to the square of the field correlation function in space 
[12, 16]. Similarly, the temporal profile of the focused 
pulse is seen in Fig. 2 to correspond to the square of 
the field correlation function in time, |F^(Ai)| 2 , where 
F° E = (E a {t')EW + Ai))/«/(f')>(J(f + At))) 1 / 2 . For 
an incident Gaussian pulse, the square of the field cor- 
relation function is equal to the intensity profile of the 
incident pulse and is independent of delay time. 

We have shown [14] previously that the contrast \i 
between the average intensity at the focus (Ip) and 
the average background intensity (h^p), in a large sin- 
gle transmission matrix of size N in steady state is, 
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/x = 1/(1/71/ — 1/N). Because the full transmission ma- 
trix is not accessible in the experiment, the density of 
measured transmission eigenvalues differs from theoreti- 
cal prediction [11,17-19]. Nonetheless, we find in steady 
state measurements [14] that when part of the trans- 
mission matrix is measured, the contrast in focusing via 
phase conjugation is given by, 

» = 1/(1/M'-1/N>). (1) 

where AI' is the eigenvalue participation number of 
the measured transmission matrix of size N' . This is a 
property of random transmission matrices and therefore 
should apply as well to transmission matrices at differ- 
ent time delays provided the field within the transmis- 
sion matrix is randomized. In Fig. 3(a), we present the 
time evolution of (M') and (fi). Near the arrival time of 
the ballistic wave, the value of M' is close to unity and 
the contrast is not described by Eq. (1). This is because 
ballistic wave is associated with the propagating waveg- 
uides modes with the highest speed and therefore the 
transmitted field is not randomized. Once the transmit- 
ted waves at the output have been multiply scattered, a 
random speckle pattern develops and the measured con- 
trast is in accord with Eq. (1). 
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Fig. 3. (Color online) Time evolution of (M') and the 
maximal focusing contrast (fi) . After the ballistic arrival, 
when [i is well described by Eq. (1), the signal is 4 orders 
greater than the noise level as indicated in the semilog 
plot of time of flight distribution shown in Fig. 3(b). 

After the arrival of ballistic waves, the value of M' is 
seen in Fig. 3(a) to increase rapidly before falling slowly. 



This reflects the distribution of lifetimes and the degree 
of correlation in the speckle patterns of quasi-normal 
modes [20-22]. Just after the ballistic pulse, transmis- 
sion is dominated by the shortest-lived modes. These 
modes are especially short lived and strongly transmit- 
ting because they are extended across the sample as a 
result of coupling between resonant centers [23-25]. Sets 
of extended modes that are close in frequency could be 
expected to have similar speckle patterns in transmis- 
sion [22], so that a number of such modes might then 
contribute to a single transmission channel. As a result, 
the number of independent eigenchannels of the trans- 
mission matrix contributing substantially to transmis- 
sion would be relatively small at early times and M' 
would be low. At late times, only the long-lived modes 
contribute appreciably to the transmission. Thus for in- 
termediate times, modes with wider distribution of life- 
times than at cither early or late times contribute to 
transmission and these modes are less strongly corre- 
lated than at early times so that M and the contrast are 
peaked. 

In conclusion, we have shown that measurements of 
spectra of microwave field transmission coefficients be- 
tween pairs of source and detector antenna positions 
make possible a detailed comparison of focusing calcula- 
tions and measurements. These results provide the con- 
ditions for optimal focusing of light and sound at a se- 
lected time. The time profile of the focused pulse is equal 
to the square of the field correlation function in time. The 
spatial contrast at the time the wave is focused depend 
only upon M' and N' of the measured transmission ma- 
trix at that time. The time evolution of maximal focusing 
is related to the range of decay times of the underlying 
quasi-normal modes and the degree of correlation of the 
modal speckle patterns. This study provides a framework 
for focusing optical pulses for bio-medical and commu- 
nications applications and for linear as well as nonlinear 
microscopies. 
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